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On post-Lie algebras, Lie-Butcher series and moving frames 
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Abstract 

Pre-Lie (or Vinberg) algebras arise from a flat and torsion free connection on a differential 



o 

Cu manifold. These algebras have been extensively studied in recent years, both from algebraic 

^^ operadic points of view and through numerous applications in numerical analysis, control 

^^ theory, stochastic differential equations and renormalization. Butcher series is an algebraic 

T-H tool used to study geometric properties of ffows on euclidean spaces, which is founded on 

Cn pre-Lie algebras. Motivated by analysis of flows on manifolds and homogeneous spaces, we 

investigate algebras arising from flat connections with constant torsion, leading to the def- 

^^ inition of post-Lie algebras, a generalization of pre-Lie algebras. Whereas pre-Lie algebras 

/^ are intimately associated with euclidean geometry, post-Lie algebras occur naturally in the 

, ^ 1 differential geometry of homogeneous spaces, and are also closely related to Cartan's method 

of moving frames. The generalized Lie-Butcher series combining Butcher series with Lie series 

are used to analyze flows on manifolds. In this paper we show that Lie-Butcher series are 

founded on post-Lie algebras. The functorial relations between post-Lie algebras and their 

M enveloping algebras, called D-algebras, is explored. Furthermore, we develop new formulas 

' ' for computations in free post-Lie algebras and D-algebras, based on recursions in a magma, 

and we show that Lie-Butcher series are related to invariants of curves described by moving 

►^ frames. 

oo 

^ 1 Introduction 



In 1857 Arthur Cayley published a remarkable paper [S] introducing the idea that differential 
operators can be described in terms of trees. This became a fundamental tool in the analysis of 
pg numerical flows through the seminal work of John Butcher |31|3j in the 1960s, where he introduced 

the Butcher group, or, in modern language, the character group of the Butcher-Connes-Kreimer 
Hopf algebra. In 1963 pre-Lie algebras appeared simultaneously from two different investigations, 
Vinberg [3D] from differential geometry (classification of homogeneous cones) and Gerstenhaber 
[T5] from algebra (Hochschild cohomology). The free pre-Lie algebra [7], described in terms of 
trees, is the algebraic foundation of B-series [TU]. B-series have evolved into algebraic tools that 
are suitable for studying a variety geometric properties of flows, such as symplecticity, preservation 
of first integrals and backward error analysis. In the late 1990s similar structures appeared in the 
renormalization theories of Connes and Kreimer [8, . Christian Brouder P] pointed out connections 
between this theory and numerical analysis. 

During the 1990s numerical integration was generalized from euclidean spaces to manifolds [21 
[24l [25l [34] . In this work it was necessary to generalize B-series to manifolds, called Lie-Butcher 
series (LB-series). Inspired by the unexpected connections between numerical analysis and renor- 
malization, the algebraic foundations of LB-series were investigated in several papers in the last 
decade [17] [T] [^ [551 HO] ■ Through this work it became clear that the foundations of LB-series 
are fundamental algebraic structures generalizing pre-Lie algebras, which may be important also 
in areas of mathematics outside of numerical analysis. 

The motivation for the investigations of this paper was to establish the algebraic foundations 
of LB-series in a similar manner to the foundation of B-series in terms of pre-Lie algebras. This 
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leads to the definition of post-Lie algebras, first found by Bruno Valette [3^ in 2007 through the 
purely operadic technique of Koszul dualisation. In this paper we show that post-Lie algebras also 
arise naturally from the differential geometry of homogeneous spaces and Klein geometries, topics 
that are closely related to Cartan's method of moving frames. Applications of moving frames in 
computational mathematics have been pioneered by Peter Olver and his co-workers |32 El [T^] • 
In this paper we will also show that post-Lie algebras and LB-series are related to moving frame 
theory and point to possible applications of moving frames in the design of numerical methods. 

2 Post-Lie algebras in differential geometry 

2.1 Algebras of connections 

In this section we motivate the abstract definition of pre-Lie, post-Lie and Lie-admissible algebras 
by considering algebras of vector fields originating from constant connections on a manifold. 

The most fundamental concept in differential geometry is connections, defining parallel trans- 
port and covariant derivations. Connections appear in various abstractions, e.g. Koszul, Ehres- 
mann and Cartan connections. To motivate pre-Lie, post-Lie and Lie-admissible algebras, it is 
sufficient to consider the simplest definition, a Koszul connection on the tangent bundle. 

Let XA4 denote the vector fields on a manifold A4 . A Koszul connection [38_ is defined as map 
V : XM X XM -^ XM such that 

^fxV = f^xV 

yx{fy)^dfix)y + fV,y, 

for any x,y € XM and scalar field /. The connection defines a (non-comniutativc and non- 
associative) R-bilinear product on XAi, henceforth written as 

x\>y -.^ Vxy- 

The torsion of the connection is a skew-symmetric tensor T : TA4 A TA4 — > TA4 defined as 

T{x,y) ^x>y -y>x-lx,yl, (1) 

where |-, •] denotes the Jacobi-Lie bracket of vector fields. The curvature tensor R: TMATA4 — > 
End{TM) is defined as 

R{x, y)z = x\> {y\> z) — y\> {x\> z) — \x,y\\> z — a{x, y, z) — a{y, x, z) -\- T{x, y) \> z, (2) 
where a{x, y, z) is the (negative) associator of the product [>, defined for any algebra {,4, [>} as 

a(x, y, z) -.^ X \> {y \> z) — [x \> y) \> z. (3) 

The relationship between torsion and curvature is given by the Bianchi identities 

&{T{T{x,y),z) + {V,T){y,z)) = &{R{x,y)z) (4) 

6((V,i?)(2/,z)+i?(r(x,y),z)) = 0, (5) 

where S denotes the sum over the three cyclic permutations of (a;, y, z). 

Torsion free connection => Lie-admissible algebra. If T = then (pjl-Q imply 

6(a(x, y, 0) -a(2/,a;, z)) = 0. (6) 

A general algebra with a product \> satisfying (rol) is called a Lie- admissible algebra. Lie-admissible 
algebras are exactly those algebras which give rise to Lie algebras by skew-symmetrization of the 
product [H], i.e. a bracket defined as 

|x,y] :=a;[>y-y [>a;, 

is a Lie bracket if and only if {A, >} is Lie-admissible. For a torsion free connection on vector 
fields, |-, •] is the Jacobi-Lie bracket. Any associative algebra is clearly Lie- admissible. A more 
general example is pre-Lie algebras. 



Flat and torsion free connection => Pre-Lie algebra. Consider a connection wliicli is botli 
flat R = and torsion free T — 0. Equation (pi) implies that 

a{x,y,z)-a{y,x,z) =0. (7) 

A general algebra with a product [> satisfying (l7| is called a pre-Lie algebra or Vinberg algebra. 
Pre-Lie algebras appear in many applications. A fundamental algebraic result is that the free 
pre-Lie algebra is the set of rooted trees with grafting as the product ^7^. This structure is the 
foundation of all classical B-series, and this was essentially known to Arthur Cayley already in 
1857 0. 

Note that a Riemannian manifold with a flat, torsion free connection is locally isometric to a 
euclidean space M" with the standard metric (3S], hence pre-Lie algebras are tightly associated 
with the differential geometry of euclidean spaces. For Lie groups and homogeneous spaces (Klein 
geometries), pre-Lie algebras are not sufflciently general to capture the basic differential geometry 
algebraically. 

Flat and constant torsion connection => Post-Lie algebra. Given a connection which is flat 
R = and has constant torsion VT = 0, then Q reduces to a Jacobi identity 6{T{T{x, y),z)) = 



and hence the torsion defines a Lie bracket (see Remark 2.2 about the negative sign). 



[x,y]:^~T{x,y). (8) 

The covariant derivation formula 'Vx{T{y,z)) — {WxT){y,z) + T{Vxy,z) + T{y,Vxz) together 
with Va;T = imply 

X \> [y, z] = [x\> y, z] + [y,x\> z]. (9) 

On the other side, (pi) together with R = imply 

[x,y]t> z:^a{x,y,z)~a{y,x,z). (10) 

This leads to the general definition of a post-Lie algebra. 

Definition 2.1. A post-Lie algebra {A, [■, •], >} is a Lie algebra {A, [■, ■]} together with a product 
\>: Ax A ^ A such that ([9|-([l0l) hold. We call [•,•] the torsion and [> the connection of the 
post-Lie algebra. 



Remark 2.2. In many applications one may naturally obtain (10) with opposite sign 

[x, y]> z^ a{y, x, z) - a{x, y, z). 

We could have defined left and right post-Lie algebras according to these sign changes. This 
terminology would make sense e.g. in the case of the Maurer-Cartan form on a Lie group. However, 
the sign in ( [lO| ) can always be absorbed into a change of the sign in the definition of the torsion 
bracket, since {yl, —[•,•]} is also a Lie algebra and (|9| is invariant under a sign change in the 
torsion. There is therefore no need to define both left- and right- versions of post-Lie algebras. 

Remark 2.3. Post-Lie algebras were introduced around 2007 by B. Vallette 02] j who found 
the structure in a purely operadic manner as the Koszul dual of a commutative trialgebra. The 
enveloping algebra of a post-Lie algebra was independently introduced about the same time under 
the name D-algebra in 29], and studied in the context of Lie-Butcher series. A basis for the 
free post-Lie algebra was presented in [27j, before the post-Lie definition was formalised. Vallette 
defines a post-Lie operad and proves that post-Lie algebras have the important algebraic property 
of being Koszul. This property is shared by many other important algebras, such as Lie algebras, 
associative algebras, commutative algebras, pre-Lie algebras, dendriform algebras etc. He also 
defines the operadic homology of post-Lie algebras and computes this for the free post-Lie algebra. 



Remark 2.4. The present differential geometric explanation of a post-Lie algebra is, as far as we 
are aware, new. Since condition ( 10 ) is expressing the flatness of the connection, while Q derives 



from the constant torsion, we initially called this structure a FCT (flat, constant torsion) algebra, 
but will henceforth adhere to the name post-Lie algebra. 

A pre-Lie algebra is a post-Lie algebra where [•, •] = 0, hence most results about post-Lie 
algebras also yield information about pre-Lie algebras. In particular, the ZJ-algebra provides a 
definition for the enveloping algebra of a pre-Lie algebra. We return to this in the sequel. 

2.1.1 Some basic results about post-Lie algebras. 

Proposition 2.5. //{A, [•,•], [>} is post-Lie, then the bracket |a;,y] defined as 

|x, y] ■.^xt>y~y>x+ [x, y] 

is a Lie bracket. 



Proof. Identifying A with Der(C/(A)), defined in Section 3.2 we get fx^yj = xoy — yox. Since o 



is associative this is a Lie bracket. D 

In the case of a flat constant torsion connection on a manifold M., the Lie bracket |-,-] is the 
Jacobi-Lie bracket of vector fields on M . 

By a modification of the product [> in A, we obtain another post-Lie algebra. 

Proposition 2.6. Let {A, [•,•], [>} be post-Lie. Define the product ► as 

x*-y = xt>y+[x,y], 

i/ien {A, —[•,•], ►} is also post-Lie. 

Proof. Since both x [> • and [x, •] are derivations on the torsion bracket, x > ■ -\- a[x, •] is also a 



derivation, for any a. A direct computation shows that (10 1 holds with a sign change, which is 



corrected by negating the torsion bracket. D 

Proposition 2.7. Let {A, [•, •], [>} be post-Lie. Define the product >- as 

x>y = x\>y-\--[x,y], 

then {A, >-} is Lie-admissible, torsion free with constant curvature 

R{x, y)z ■.= a^{x,y,z)- a^ {y, x, z) = - - [[x, y],z]. 

Proof. Lie-admissible follows from x >- y — y >- x = la;, j/]. The curvature follows from a lengthy 
but straightforward computation. D 

Remark 2.8. In the case of vector fields on a Lie group, discussed below, [> and ► come from 
the right and left Maurer-Cartan form and >- is the Levi-Civita connection. 

2.2 Lie groups, homogeneous spaces and moving frames 

A more general view on torsion and curvature appears in the theory of G-structures and g-valued 
forms on a manifold. This is the foundation for Cartan's method of moving frames, which Pe- 
ter Olver and his co-workers have developed into a powerful tool in applied and computational 
mathematics [33l ES . 



2.2.1 Post-Lie algebras and numerical Lie group integrators 

We will briefly review the fundamentals of numerical Lie group integration, as formulated in 
Let G be a Lie group with Lie algebra g, and let 

X: G X M ^ M, {g,u)\^g-u 

be a transitive left action of G on a homogeneous space AA , with infinitesimal generator 

d 



K:qxM^TM, {v,u)^ 



exp{tv)-ueTuM. (11) 

"'- t=o 

Let ri'^(A^, g) be the space of g-valued k-forms on M, in particular 51° (A^, g) is the space of maps 
from Ai to g. Any x £ rf{M,g) generates a vector field X E XAi as 

X{u)=Kixiu),u), (12) 

which by abuse of notation is written compactly as X = X-t{x), where A* : 51°(A^,g) — > XA4. 

Remark 2.9. Most Lie group integrators for the differential equation u'{t) = F{u(t)), where 
u{t) G M. and F g XM, are based on rewriting the equation as u'(t) = X^,f{u{t)) for / G 
il'^{A4,g). It is important to note that if the action of G is transitive but not free on Ai, then 
A* : il*'(A^,g) — ?► XAi is surjective but not injective. The freedom in choice of a / to represent F 
is given by the isotropy (stabiliser) subgroup of G at a point u E A4. Different choice of isotropy 
can lead to significantly different numerical integrators. As pointed out by Lewis and Olvcr ^fS] . 
moving frames is an important tool in the study of isotropy choice for Lie group integrators. We 
return to this in the sequel. 

Numerical analysis of Lie group integrators is intimately related to post-Lie algebras, due to 
the following result. 

Proposition 2.10. Let Ad he acted upon from left by a Lie group G with Lie algebra {g, [•, -Jj,}. Let 
the Lie bracket [■,■]: n°[M,g)y.n°{M,g) -)■ n°{M,g) and the product > : n°{M,g)xn'^(M,g) -^ 
ilP {A4 , g) be defined pointwise at u E M as 

[x,y]{u) = [x{u),y{u)]^ (13) 

x[>y — X^{x){y) (the Lie derivative of y along X^,{x)), (14) 

then {fl'^{A4,g), [•, •], [>} is a post-Lie algebra. 

Proof. This can be verified by a coordinate computation. Let {cj} be a basis for g and dj = A*(ej) 
the corresponding right invariant vector fields on A4. Note that A*([ej,efe]) — —ldj,dkl, where 
the right hand side is the Jacobi-Lie bracket of vector fields. Letting x{p) — J^j ^''ip)^j ^^^ 
y{p) — Sfc y''{p)^ki where x^ and y'^ are scalar functions on A4, we obtain 

[x,y] = '^^'y''[ej,ek] 

j,k 

x\>y = ^x^dj{y'')ek. 
j,k 

The post-Lie conditions follow by a straightforward computation. D 

A similar result, formulated in the enveloping algebra, can be found in [291 Lemma 3]. 

The connection [> leads to a parallel transport of vector fields (f>*x{p) = x{4>{p)) for x E 
il^{A4,g) and (p G Diff(A^). This parallel transport is used in numerical Lie group integrators to 
collect tangent vectors to a common base point in order to compute the timestep of the algorithm. 
The parallel transport is clearly independent of paths, since the connection is flat, and it is given 
algebraically as the exponential of >. 



2.2.2 The Maurer-Cartan form 



The (left) Maurer-Cartan (MC) form on a Lie group G is a g-valued one- form u E 51^(G,g), 
defined as the map w: TG — )• g moving v G TgG to g by left translation 



uj(v) 



TLg-iv, 



where Lgh — gh is left multiplication in the grouprl 

The Maurer-Cartan form defines a linear isomorphism utp : TpG — )• g and hence defines an 
isomorphism between f2°(G,g) and vector fields XG. For the right Maurer-Cartan form, this 
isomorphism is given by A* defined in (ff |-(f2), when M = G. For the left MC form, the 



corresponding isomorphism p.^ : fl'^{G,g) — > XG is given as the infinitesimal right action 



p*ix){g) 



d_ 
di 



/exp(te(5)). 



t=o 



The Maurer-Cartan form satisfies the structural equation 



duj + -uj A uj 



0. 



(15) 



(16) 



On a general (connected, smooth) manifold Ai , the existence of a g-valued one form which is an 
isomorphism on the fibre and satisfies ( 16 ) implies that Ai can be given the structure of a Lie 



group (up to a covering) [36l Theorem §8.8.7]. Thus the Maurer-Cartan form is fundamental in a 
differential geometric characterization of Lie groups. 



The curvature of w G ^^ (G, g) is given as R = duj + ^lu A oj E fl^ (G, g) , and ( 16 1 is therefore 
a flatness condition R = 0. Taking 9 = —uj as a solder form, we compute the torsion form 
Q^de + LuAd 



-^cjAojE n^{G, g). This yields 

e(x,r) = -Kx),^(y)]g. 

Therefore, the Maurer-Cartan form is flat with constant torsion. 

Proposition 2.11. Given a Lie group G and the inverse of the (left) MC form p* : fl'^{G, g) 
in (15), then {il*'(G, g), —[•,•], ►} is a post-Lie algebra, where 



XG 



x*-y := p*{x){y) 
[x,y]{g) := [x{g),y{g)]Q. 



The product x>- y is related to the connection of the right MC form x\> y, see {14h o,s 

X ^ y = X \> y + [x,y]. 



(17) 



Proof. From the right MC form we get the post-Lie algebra {fl'^{G, g), [•,•], [>} (Proposition 2.10 ). 
The correspondence between a right trivialized x and the corresponting left trivialized x in G 
i^°{G,g) is given as x {g) = uj{X^,x){g) = Ad gx{g), from which (17) follows by differentiation. 
Hence, by Proposition 2.6 {ri'^(G,g), — [•, •],►} is post-Lie. D 



Remark 2.12. Since p* : i7°(G, g) — ?► XG is an isomorphism, we can equivalently define ► : XGx 
XG — > XG as X ► y = p* {X{lu{Y))). This is a flat Koszul connection on XG with torsion 
[X, F] — —p^,[uj{X),uj{Y)]g. The vector fields on a Lie group with this connection and torsion is 
a prime example of a post-Lie algebra. 



^The MC form can also be defined by right translation, but the left form is more convenient for moving frames. 



2.2.3 Homogeneous spaces, a leftist view 

We recall some aspects of the Klein-Cartan geometry of a homogeneous space M from [35] . Given 
a transitive left action ■: Gx M. ^ Jvl and an arbitrary point o e A^. Let H = {h G G \ h-o ~ o} 
be the isotropy subgroup at o with Lie algebra () < g. Define 

G x_H := (G X fl) / ~ , 

where {g,v) ^ {gh, Ad h-iv) for all h E H. Define the map 

Pm- G XhQ^ M, {g,v)^ gexp{v) -o, (18) 

and its differential with respect to the second variable pM* : G Xh g -^ TAi as 



PM*(g,v) := g-^ 



gexp{tv)-oeTg.oM. (19) 



Since pM^,{g,v + v-^) ~ PM^{g,v) for all v^ G f), it follows that pM* induces a smooth quotient 
mapping 

m*- GxHQ/i)^TM,{g,v + b,)^ pM^{g,v), (20) 

where g/f) denotes the quotient as vector spaces. In [35] [Prop. 5.1] it is shown: 

Proposition 2.13. 7?/v7* defines an isomorphism G Xu g/f) ~ TM as vector bundles over A4. 

Thus, tangents in Tg.oM are uniquely represented as {g,v) G G Xh fl/f), while finite motions 
are not invariant under change of isotropy; in general pM{g,v) ^ PM{g,v + w^) for ti-'- G f). In 
order to fix a choice of isotropy, it is useful to discuss the notion of gauges. This is essentially the 
same concept as gauges in field theories of theoretical physics. 

A Cartan gauge is a local section of the principal iJ-bundle tt: G — > A^,.g f-> 5-0, i.e. a map 
a: U C Ai —?' G such that noa = Id on an open set U. We denote this (ct, U). 

The Darboux derivative of a map / : A^ — >■ G, denoted w/ : TM — > g, is defined as the puUback 
of the Maurer-Cartan form w on G along /, 

^f -^ f*^ = wo/*, (21) 

where /* : TAi —J- TG is the differential of /. This will always satisfy the Cartan condition 

ujf + -ujfAujf^O, (22) 

and thus it is a flat g-valued 1-form ujf G i^^{A4,g). We call / the primitive oiujf. A generalization 
of the fundamental theorem of calculus [36j states that, locally, a one form 6 G i^^{Ai,g) has a 



primitive /: A^ — > G if and only if 6 is flat, i.e. iff 9 satisfies (22). Furthermore, the primitive is 
determined uniquely up to left multiplication by a constant c (z G: if Uf — uj: then f — cf for 
some constant of integration c G G. 

The Darboux derivative of the Cartan gauge uia- G fl^{M,g) is called an infinitesimal Cartan 
gauge. Since a;^ is flat, we can uniquely recover a from iOa- (the integration constant is found by 
the condition that ct is a section). 

Proposition 2.14. A Cartan gauge {a, U) with derivative Ua-, defines a map cr* : TA4 — >■ G x^ g 

as 

a^{v) = a{Trv)xHUJa{v), (23) 

which extends to (f* : TAi — )• G x^ g/f) by composition with the projection g/t) 1—^ f) 

a,:TM^GxHg^GxHg/i}. (24) 

The maps a 1, andpM^ are locally inverse vector bundle isomorphisms, (7*0^0^7* = Id andpM^.^o'* = 
Id on their domains of definition. 



Proof. This is a straightforward computation. D 

Infinitesimally all gauges are equivalent, but how does a change of gauge affect the finite 
motions induced by pM? A retraction map is a smooth, locally defined map TZ: TAi — > A^ such 
that: 

• 7?.(ti) = TTi\iv if and only if i; is a 0-tangent, where ttm ■ TAi — ?► Ai. 

• TZ'{0) — Id (the identity on the tangent fibre). 

Retractions provide a useful way of formulating numerical integration schemes f6l . For any Cartan 
gauge (cr, U) there corresponds a retraction map TZa ■ TU — > Ai defined as 

TZa{v) — pM°cr*{v) = a{TTv)exp{uJc{v))-o, iorvGTU. (25) 

If {a, U) and (cr', U) are two gauges on [/ C A^, then (j'{u) — a{u)h{u) for some smooth h: U ^ H. 
The corresponding infinitesimal gauges w^ and uta' are related as [36] [p. f68] 

J^^ kAh-iuj„+ujh, (26) 

where LOh G ^^{U, f)) is the Darboux derivative of h. Thus, the modified retraction is given as 

TZa' = cr(7rw)exp(a;cr(w)+t^/i(u))-o. (27) 

If h is constant then ujh — ^ and the two retractions are equal, but generally they differ. For 
numerical integration it is important to choose the gauge in a good manner. We return to this in 
the discussion of moving frames below. 

Let r(G y-H q) denote sections of the bundle G x^ g ^- A^. This space can be identified with 
the subspace r2^(G,g) C n^{G,Q) defined as 

n°j{G,g) :={xe n"{G,g) \ x{gh) = Ad„-^x{g) }. (28) 

Proposition 2.15. There is a 1-1 correspondence between x £ il'jj{G,g) and X G r{G x^ g), 
defined locally by a Cartan gauge (cr, U) as 

X{u) = a{u) Xh x{a{u)), for u G U C M. 

The map x i-^ X is independent of the choice of section a. 

Proof, x I— >■ X is independent of a since a{u)h x x{a{u)h) ~ cr{u) x Ad hx{cr{u)h) — cr{u) x x{a{u)). 
Given X we recover x on the section a{U) C G, and reconstruct x on the fibres of if — >■ G ^- A^ 
^y x{gh) = Ad h-ix{g). D 

Thus, to sum up, we have an identification il'jj{G,Q) ~ r(G x^ g). An X £ r(G x^ q) 
corresponds to a vector field pm*oX £ XAi and a diffeomorphism pmoX £ Diff(A^). Furthermore, 
given a Cartan gauge (cr, U), we can map a vector field Y £ A'A^ to cr*y £ r(G Xh g), but this 
map depends on the gauge a. Finally, there is a post-Lie algebra also associated with this view of 
homogeneous spaces: 



Proposition 2.16. {f75f(G,g), —[•,•], ►}, with ► and [■,■] defined in Prop. 2.11, is post- Lie 



Proof. A straightforward computation reveals that (x ► y){gh) = Ad^-i(a; ► y){g) and 
[x{gh),y{gh)] = Ad fi-i[x{g),y{g)]. Hence ► and [•, •] are well defined on the subspace Q°jj{G,q). 

D 

What is the parallel transport in this post-Lie algebra? Note that x £ ft^{G,g) defines a 
vector field p^,x £ XG with a fiow $f : G —?' G satisfying 

^t{gh) = <^tig)h, for all h£H. (29) 

For such a diffeomorphism we define $* : il'jj{G, g) -^ il''^{G, g) as 

<^*y{g) := y($(g)). (30) 



Proposition 2.17. For x,y E fl'fj{G,g) define the exponential series 

exp{x>-)y := y + X >■ y + -{x >■ {x >■ yj) + -{x >■ {x >■ {x >■ y))) H . (31) 

2 D 

Parallel transport in il^(G, g) is given as 

exp{x>-)y = $*y, 

where ^* : G ^ G is the t = 1 flow of p*x G XG. 

Remark 2.18. The post-Lie structure captures algebraically both infinitesimal aspects of homo- 
geneous spaces and also finite motions such as flows and parallel transport. It is therefore clear 
that the post-Lie structure cannot carry over to the quotient space F (G Xh fl/f)) — XM, both ► 
and [•,•] change under an isotropy change x i-^ x + x-^ where XjX-^ E n'jj{G,Q) and x'^{g) S (]. 
Hence this post-Lie structure does not carry over to XA4, except when the action is free. 

Remark 2.19. In Section [2.2.1[ we described a post-Lie algebra on a homogeneous space A4 
derived from the right MC form. This has been the basis for almost all work on numerical 
integration schemes on homogeneous spaces since |30j . In the present section, we have detailed 
an alternative post-Lie structure on A4 derived from the left MC form. We believe that this 
formulation should lead to better geometric integration algorithms on symmetric spaces. Also, it 
seems as this formulation is better suited for combination of numerical Lie group integration with 
moving frame algorithms. 

2.2.4 Moving frames 

Moving frames provide an important tool for choosing gauges that are naturally derived from the 
geometry of e.g. a differential equation or other geometric objects such as curves and surfaces 
in a homogeneous space (Klein geometry). Let A^ be a manifold acted upon from left by a Lie 
group G with Lie algebra q. We do not require the action to be transitive, so A4 needs not be a 
homogeneous space. 

Definition 2.20. A left moving frame is a map a: Ai -^ G such that 

a{g ■ u) — ga-{u) for all g e G and u G A^, 

a right moving frame is a map r:M^G such that 

r{g ■ u) — r{u)g^^ for all <? € G and u £ M.. 

If r is a right moving frame then a (it) = r(u)~^ (inverse in G) is a left moving frame. Moving 
frames exist if and only if the G action on A^ is free and regular. In that case, moving frames can 
be constructed (locally) as follows [55] : 

1. Choose a submanifold /C C A^ which is transverse to the G orbits and of the maximal 
dimension p — dim(il/) — dim(G). Locally, there is one point in /C for each orbit, and each 
orbit intersect /C in one point. In coordinates, K. is often chosen by setting d — dim(G) of 
the coordinates to constant values. 

2. A right moving frame r is found by solving the normalisation equations r{u)u e /C for r{u). 

3. A left moving frame is obtained by inverting r. 

If the action is not free, there is a standard procedure of obtaining a free and regular action by 
prolongation of the group action, i.e. we extend A/J to a jet-space J'^(A^), which is the geometrical 
way of saying that we consider the space of all curves in M. represented by Taylor expansions up 
to order k. The prolongation of the group action is the natural induced action of G on (Taylor 
expansions) of curves. Coordinates on the jet-space are given by the (higher order) derivatives of 



curves. We can always obtain a free regular action (and thus a moving frame) by prolongation. 
Thus, by this construction we find a left moving frame a: J^{M.) — ?► G. 

Moving frames are closely related to Cartan gauges on a homogeneous space M.. Let G act 
transitively on AA with an isotropy subgroup i/, and let ct: A^ — ?► G be a local section of the 
bundle it: G ^f M.. Note that cr is a section if and only if 

a{g-u) = ga{u)h{u) for h: M^ H. (32) 

Thus it is a left moving frame, up to isotropy. Such a map is also called a partial moving frame |17] . 
Thus, the theory of moving frames (full and partial) provides geometric ways of construct- 
ing sections {a^ C/) of G — > A^ , hence also geometric ways of fixing isotropy through the map 



a^ : TM -^ T{G Xh q) ~ fl%{G,Q) in (23 1. On n%{G,g) we have aU the tools we need to do nu- 
merical integration and analysis of numerical integration schemes. The details of such algorithms 
is subject to future research. We see at least two useful ways to proceed in choosing a. 

• By prolongation of the group action we can obtain a full moving frame a: J^{M.) — > G. To 
solve a differential equation u' — F{u), F £ XAi, we must prolong also F to the jet-bundle. 
This should be a very attractive numerical method in cases where we can compute the fc-th 
derivatives of F, either by computer algebraic means, or by automatic differentiation. 

• In the case where AJ is a symmetric space, there is a canonical choice of section ct : A^ — >■ G. 
In this case there exists a canonical splitting q = t) ®i, where f) is a sub algebra and t is a 
Lie triple system (LTS). The infinitesimal gauge uj^ takes values in i, and thus exponentials 
need only be computed on the LTS. Efficient algorithms for computing exponentials on an 
LTS are discussed in |41j . This theory opens up the possibility of new classes of numerical 
integration on symmetric spaces. 

3 The algebraic structure of post-Lie and D-algebras 

In this chapter we discuss the algebraic structure of general post-Lie algebras {A, [•, ■],[>}. Various 
aspects of this theory can also be found in |571 [321 HH [S] ■ Here we develop the core theory from 
the axiomatic definition of a post-Lie algebra, and establish the functorial relationship between 
post-Lie algebras and the enveloping _D-algebra. Moreover, we find that a magmatic view of planar 
trees and forests give rise to recursive formulas for various algebraic operations, which simplify 
computer implementations. 

3.1 Free post-Lie algebras 

In [7] Chapoton and Livernet gave an explicit description of the free pre-Lie algebra in terms of 
decorated rooted trees and grafting. In this section we will see that there is a similar description 
of the free post-Lie algebra. In fact, we will show that the free post-Lie algebra can be described 
as the free Lie algebra over planar rooted trees, extended with a connection given by left grafting 
of trees. Furthermore, we will relate post-Lie algebras to D-algebras, studied in connection with 
numerical Lie group integration ^\ [52] . The universal enveloping algebra of a post-Lie algebra is 
a D-algebra, and the post-Lie algebra is recovered as the derivations in the D-algebra. 

Trees. Let C be a set, henceforth called colors. We define Tc the set of all planar (or orderedlj 
rooted trees with nodes colored by C. Formally we define this as the free magma 

Tc :=Magma(C). 

Recall that a magma is a set with a binary operation * without any algebraic relations imposed. 
The free magma over C consists of all possible ways to parenthesize binary operations on C. We 

^Trees with different orderings of the branches are considered different, as when pictured in the piane. 
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identify Magina(C) with planar trees, where the nodes are decorated with colors from C. On trees 
we interpret • as the Butcher-product [1] : ri * r2 = r is a tree where the root of the tree ti is 
attached on the left part of the root of the tree T2. For example: 



jJ.^ 



= (•*o) •((••(•••)) •( 



If C = {•} has only one element, we write T :— T{,}. The first few elements of T are: 



.,:Yi,^,V,^/,Y, 



Note that any r G Tc has a unique maximal right factorization 

T = Ti * (t2 * (• • • (r^ * c))), where c £C and ri , . . . , rj. e Tg . 

Here c is the root, k is the fertility of the root and ri, . . . , r^ are the branches of the root. Let k 
be a field of characteristic zero and write k{Tc} for the free k- vector space over the set Tc, i.e. all 
k-linear combinations of trees. We define left graftin^\> : T^xTc — > k{Tc} by the recursion 

T [> c :— T -k c, for c € C 

T [> (ti * (t2 • (• • • (Tfe :*r c)))) := T * (ti • (t2 * (• • • (rfc • c)))) 

+ (r[>Tl)*(T2*(---(Tfe*c))) 

+ ri*((T[>r2)*(---(Tfe*c))) 



+ Tl • (t2 • (• • • ((r \> Tfc) • c))). 

Thus Ti [> r2 is the sum of all the trees resulting from attaching the root of ti from the left to all 
the nodes of the tree T2. Example: 



J>Y: 



V+Y+Y. 



Free Lie algebras of trees. Let g — Lie(Tc) denote the free Lie algebra over the set Tc 
For C = {•}, a Lyndon basis is given up to order four as [27]: 



Lie(Tc) = k 



.,I,I,YI.],I,Y,V[L1,\^, 



^JY.],[I. 



7 ^ 5 ^: ^ 1 ^ ; 



Proposition 3.1. Let the free Lie algebra g = Lie(Tc) be equipped with a product > : x g — )• g, 
extended from the left grafting defined on Tc in ^3!^ as 

(34) 
(35) 



u\>[v,w] — [u\> v,w] + [v,u\> w] 
[u,v]\>w — a{u,v,w) — a{v,u,w) 



for all u,v,w (li g. Then {Lie(Tc), [■,■], f>} is post-Lie. 

Proof. Since any u,v,w G g can be written as a sum of trees and commutators of trees, the 
connection is well-defined on g. By construction it satisfies the axioms of a post-Lie algebra. D 



^Various notations for similar grafting products arc found in the literature, e.g. uS> v = u[v] 
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Free post-Lie algebras. Proposition |3.1| shows that the free Lie algebra of ordered trees has 
naturally the structure of an post-Lie algebra postLie(C) := {Lie(Tc), [•, •], >}. We call this the 
free post-Lie algebra over the set C for the following reason: 

Theorem 3.2. For any post-Lie algebra {A, [•, •],[>} and any function f : C ^ A, there exists a 
unique morphism of post-Lie algebras T : postLie(C) — > A such that T{c) = f{c) for all c ^ C. 

Proof. We construct J- in two stages. First we show, using [>, that / extends uniquely to a 
function Tt^ '■ T^ — > ^. Then by universality of the free Lie algebra, there is a unique Lie algebra 
homomorphism T: Lie(Tc) — > A. We show that this is also a homomorphism for the connection 
product \>. 

To construct the extension to Tc we first observe that the magmatic product t -k t' on Tc 
(the Butcher product of two trees) can be expressed in terms of left grafting \> . This is done by 
induction in the fertility of r'. For fertility 0, i.e. r' = c G C, we have T-k c = t\> c. For fertility k 
we write t' — Ti-k (t2 * (• • • (r^, 7k- c))) and find from p3[) 



T • r' = T [> r' - (r [> Ti) * (t2 * (• • • (r^ • c))) (n * (t2 * (• • • (r [> r^ * c))). 

In the right hand side of the equation, the fertility of any term on the right hand side of a ^-product 
is smaller than k, which completes the induction. The fact that Tc is freely generated from C by 
the product * ensures that J-t^ is uniquely defined by 

J"tc(c) = /(c) for all cGC 

and hence that also J^: Lie(Tc) —>■ A is uniquely defined as a Lie algebra homomorphism. 

Finally, by induction on the length of iterated commutators, we see that J-{ut>v) = J-{u)\>J^{v) 
for all u,w e Lie(Tc): li u,v G Tc this holds by construction. Assuming that J^(u l> v) — 



J^{u)\>J-(v) whenever u and v are iterated commutators of length at most fc, we find by using ( 34 ) 



(351 that J"([M,ri][> [w,T2]) = J"([u, n]) [> J"([u, tj]) for allTi,r2 € Tc. D 



Proposition 3.3. Let postLie(C) be graded with the number n counting the number of nodes in 
the trees. Then 



1 -^r^ ,n^ (2d 



(iim(postLie(C)„) = -— >, m( 



2n ^-^ d \ d 



\c\ 



where fi is the Mobius function. For \C\ = 1 the dimensions are 1, 1, 3, 8, 25, 75, 245, .... See also 
J|3 A022553]. 



Proof See [1S| and [H]. D 

Remark 3.4. The same dimensions also appear for the primitive Lie algebra of the Hopf algebra 
CQSym (Catalan Quasi-Symmetric functions) (31j. We believe that also this theory is founded on 
post-Lie algebras. 

3.2 Universal enveloping algebras 

D-algebras. In Section |4] we describe certain algebraic structures that occur naturally in the 
study of numerical integration methods on manifolds |29j . Central to this work are algebras of 
derivations, called Z?-algebras. We will see that post-Lie algebras relate to D-algebras similarly to 
the relationship between a Lie algebra and its universal enveloping algebra. 

Definition 3.5 (D-algebra ,29J). Let B be a unital associative algebra with product u, w i— ^ uv, 
unit I and equipped with a non-associative product • [> • : i? (g) B — >■ _B such that I\> v = v ioi all 
V £ B. Write Der(_B) for the set of all u E B such that u [> • is a derivation: 

Der(i3) — {u £ B \ ul> {vw) = {ul> v)w + v{u > w) for all v,w £ B}. 
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B is called a D-algebra if the product u,v i-^ uv generates B from {I,Der(i3)} and furthermore 
for any u G Der(_B) and any v,w £ B we have 

vl>u e Der(B) (36) 

(uv) \> w — u\> (v \> w) ^ (u\> v) \> w. (37) 

Proposition 3.6. If B is a D-algebra then the derivations Der(i?) form a post-Lie algebra, with 

torsion [u, v] = uv — vu and connection [>. 

Proof. If u, V e Der(i3) we note that 

(uv — vu) \> ■ = u\> {v\> ■) — v\> {u\> ■) + {u\> v)\> ■ — {v\> u)\> ■. 

The first two terms on the right is a commutator of two derivations and is therefore a derivation. 
The last two terms are derivations separately. Hence, [u,v] £ Der(i?) and {Der(i3), [•, •]} is a Lie 
algebra. The other axioms of being post-Lie follows easily from the definition of a D-algebra. D 

Universal enveloping algebras. Let {A^ [•,•], [>} be an post-Lie algebra, and let U{A) be the 
universal enveloping algebra of the Lie algebra {A, [•,•]}. By the Poincare-Birkhoff-Witt theorem 
we can embed A as a linear subspace of U{A), such that [m, v] — uv — vu. The embedding is also 
denoted by A. The product \> on A can be extended to U{A) according to: 

I[>u = V (38) 

u \> (vw) — {u \> v)w + v{u \> w) (39) 

{uv) \> w = u\> {v\> w) — {u\> v)\> w, (40) 

for all u G ^ and v,w E U{A). 



Proposition 3.7. Equations !(38\l-{40) define a unique extension of\> from A to U{A). With the 



non-associative product [>, U{A) is a D-algehra with derivations Dcv{U {A)) = A. 

Proof. See (161 Theorem V.l] for a proof that a derivation on a Lie algebra A extends uniquely 
to a derivation on U{A). This justifies the extension on the right ([39|. The extension on the 



left, given by (38) and (|40p, is compatible with the the embedding [u,w] ^-^ uv — vu due to the 



flatness condition (10) for post-Lie algebras. From the PBW basis on U{A) it follows that these 



equations extend > uniquely to all of U{A) also on the left. It is clear that A C Der(C/(^)). 



To check that A — Der([/(yl)) we verify from (p?8|-(40) that I is not a derivation and that 



Mi,W2 G Der(C/(^)) => ui?i2 <^ Der([/(v4)), thus Der(C/(A)) cannot be larger than A. D 

Definition 3.8 (Universal enveloping algebras). We call U{A) equipped with this D-algebra 
structure \> the universal enveloping algebra of the post-Lie algebra A. 

Proposition 3.9. For any D-algebra B and any post-Lie morphism f : A ^f Der(_B) there exists 
a unique D-algebra morphism T : L!{A) — > B such that J-{u) — f{u) for all u £ A. 

Proof. J- is uniquely defined as a unital associative algebra morphism. It remains to verify that 
T(u [>«)=: .7^(u) > J'iyi)- f^(^) has a grading by the length of the monomial basis of PBW. 



Using (38|-(40) it follows by induction in the grading that T(u \> v) — T(u) \> J-(v). D 



Remark 3.10. The preceding results establishes a pair of adjoint functors between the categories 
of D-algebras and post-Lie algebras: 

[/(•) : post-Lie : Z^ D-alg : Der(-). 

In other words, there is a natural isomorphism 

HompostLic(Dcr(A),B) ^ HomD(A, C/(B)). 
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Free D-algebras. A direct consequence of Theorem |3.2| and Proposition |3.9| is the foUowing 
characterization of a free D-algebra: 

Corollary 3.11 ([5Sl Proposition 1]). The algebra Dq '■— C/(postLie(C)) is the free D-algebra 
over the set C, i.e. for any D-algebra B and any function f : C ^>- Der(i3) there exists a unique 
D-algebra morphism J- : Dc — > D such that J-{c) — f(c) for all c (z C. 

The unital associative algebra of Dc is C/(Lie(Tc)), which by the Cartier-Milner-Moore theo- 
rem is the free associative algebra over Tc- I.e. it is the noncommutative polynomials over rooted 
trees: Dc = k(Tc) = k{Fc}, where k{Fc} denotes the free vector space over the set of ordered 
forests. Fc := Tj consist of all words of finite length over the alphabet Tc, including the empty 
word I. For C = {•} these are 

F = J I, •, •«, I •••, I., .1 Y, 



We can create a tree from a forest u! by applying the operator B J" : Fc — ;► Tc , attaching the trees in 
uj onto a common root labelled by c € C and we can create a forest from a tree using the operator 
B~ : Tc — 7> Fc removing the root. The concatenation product wi,W2 '~^ (^11^2 is the associative 
operation of sticking shorter words together to create longer words. 

Summarizing, the free D-algebra Dc is the vector space of forests k{Fc} with unit I, concate- 
nation product and the left grafting product \> defined on trees in (p3|) and extended to forests 



by (38)-(40|. This free D-algebra carries a Hopf algebra structure, closely related to the Connes- 



Kreimer Hopf algebra, to be discussed in the sequel. 

The composition product o on D-algebras. A dipterous algebra |TS] is a triple {i3,o,[>}, 
where _B is a vector space and o and \> are two binary operations on B satisfying: 

xo[yoz) = {xoy)oz (41) 

x> {y> z) — (xoy) \> z (42) 

for all x,y,z G B. Let i? be a D-algebra with concatenation x^y t-^ xy and connection product 
x\> y. Define a product o : _B x i? — >• i? as 

lay = y 

xoy :— xy + X \> y (43) 

{xy)oz :— xo[yoz) — (x \> y)oz for all x e Der(i?), y,z ^ B. 

Proposition 3.12. If B is a D-algebra then {B, o,[>} is a dipterous algebra. 

Proof. Proof by induction in the grading on B provided by the PBW basis. D 

The product x,y t-^ xoy will be referred to as the composition product, while x,y 1-^ xy is called 
either concatenation of frozen composition, due to the interpretation for differential operators 



on manifolds. Let A = n^{M,g) be the post-Lie algebra defined in Proposition 2.10 and let 
B = U{A) — J7"(A^, t/(g)). For f,gGB the frozen composition is {fg){p) = f{p)g{p), where we 
'freeze' the value of / and g in a point p (z A4 and obtain the product from U{q). The composition 
f,g^-^ fog, on the other hand, corresponds to the fundamental operation of composing two 
differential operators on M.. For f,g^ Der(i3) we have fog = fg + f\>g, splitting the composition 
in a term fg where g is 'frozen' (constant) and a term f \> g where the variation of g along / is 
taken into account. 

On the free D-algebra Dq the composition is computed on two forests uJi,u!2 G Fc as ( \'29\ 
Definition 2): 

LLI1OUJ2 ^ B^{uji\> B^{uj2)). (44) 

We call this the planar Grossman-Larson product, since it is a planar forest analogue of the 
Grossman-Larson product [15] of unordered trees appearing in the Connes-Kreimer Hopf algebra. 
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3.3 Hopf algebras 

Hopf algebraic structures related to the free D-algebra Dc = C/(postLic(C)) has been studied 
in [211 mini]. These Hopf algebras can both be seen as generalizations of the shuffle-concatenation 
Hopf algebras of free Lie algebras as well as of the Connes-Kreimer Hopf algebra, which is closely 
related to pre-Lie algebras [7]. 

Shuffle product. From the classical theory of free Lie algebras, it follows that the deriva- 
tions Der(Dc) can be characterized in terms of shuffle products. Define the shuffle product 
LU : Dc (8) Dc ^- Dc on the free D-algebra Dc by iLUw = uj = wLUl and 

(riWi)Lij(r2a;2) = ti{uji{1AT2UJ2) + r2(Tia;iLiJW2) 

for Ti,T2 € T, uji,(jJ2 € F. Let (•, •) be an inner product on Dc defined such that the forests form 
an orthonormal basis, and let the coproduct Alu : Dc — >■ Dc Dc be the adjoint of LU. 

Proposition 3.13. The free D-algebra Dc has the structure of a cocommutative Hopf algebra 
T-L'f^ = {k{Fc}, e, o, ry, Alu, 5'} with product being the planar Grossman-Larson product o defined 
in l44]li ihe coproduct Alu is the adjoint of the shuffle and the unit rj and counit e are given as 



7/(1) = I 

e(I) = 1, e(a;)=0 for all u: (^¥c\{l}. 
The primitive elements are PrimiT-L'j^) — Der(Dc). The antipode S is defined in \29^ . 

Proof. The Hopf algebraic structure (for the dual of "H^) is proven in [29j. Characterization of 
the primitive elements follows from the free Lie algebra structure |35) . D 

The Hopf algebra "Hat, a magmatic view In the study of numerical integration on man- 
ifolds it is important to characterize flows and parallel transport on manifolds with connec- 
tions algebraically. It is convenient to base this on the dual Hopf algebra of Ti-'j^- Let T-Lm = 
{k{Fc}, e, LIJ,?7, Ao, S*} be the commutative Hopf algebra of planar forests, where the product is 
the shuffle product LU and the coproduct Ao the adjoint of the planar Grossman-Larson product. 
Various expressions for Ao and the antipode S are derived in [29]. Our definition of Fc and T-Lm is 
rather involved, going via trees and enveloping algebras extending [> from derivations, introducing 
the dipterous composition o and dualizing to obtain Aq. However, both Fc and the Hopf algebra 
■Htv can alternatively be defined in a compact, recursive manner. We will review this definition, 
which is the foundation for a computer implementation of T-Li^ currently under construction. 

Definition 3.14 (Magmatic definition of Fc). Given a set C we let {y.c}cec be a collection of 
magmatic products, without any defining relations. Let I denote the unity and we define Fc as 
the free magma generated from I by the magmatic products. 

This definition is related to our previous definition of Fc by interpreting tui x ^ ci;2 in terms of 
forests as 

uji -x c UJ2 ^ ujiB+ {u:2) (45) 

for all 101,102 € Fc, c <E C. Thus, e.g. for c = o we have I x^, I = o, and 



• X r !• = • B . 



Any w G Fc \{I} can be written uniquely as w = ulXcUJh, where c € C is the root of the rightmost 
tree in the forest. We call u>l and cor the left and right parts of to and c the right root. 

Definition 3.15 (Shufile product). The shuffle product LU: k{Fc}«)k{Fc} -S- k{Fc} is defined 
by k-linearity and the recursion 

IlllLo = Loiul = Lo, for aU w e Fc, 

(46) 
vlULo = {vLlUio) XcVr + {vIUlol) XdiOR, for V = vl x^vb., uj = lol x^lor. 
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Definition 3.16 (Coproduct.). The coproduct Ao : k{Fc} -^ k{Fc}(8)k{Fc} is defined by k- 
linearity and the recursion 

Ao(a;) = w(g)I + Ao(a;i)LiJx<iAo(w_R,), for lj = u^ Xdi^R, 
where LLlx^ is the shuffle product on the left and the magmatic product x^ on the right: 

(ui(g)U2)LiJXrf(lii®U2) := {ui\±iVi)(^{u2 XrfW2). 

Proposition 3.17 ([IH])- T^n — {k{Fc}, e, LU, 77, Aq, 5} is a commutative Hopf algebra. 
The Hopf algebra Hn is the setting for Lie-Butcher series. 

4 Lie— Butcher series and moving frames 

Lie-Butcher series are formal power series for flows and vector fields on manifolds, which combine 
Butcher series with Lie series. An extensive overview of this field can be found in [52]. We will 
briefly review the basic definitions and point to some relations between LB-series and moving 
frames, which have not been discussed elsewhere. 

Definition 4.1 (Lie-Butcher series). Let H*j^ = Homk(HAf,k) denote the linear dual space of 
Hn- An element a G Wjq is called a Lie-Butcher series. We identify a with an infinite series 

via a dual pairing (•, •) : H^ x "H^r — > k defined such that 

a{uj) — (a,w) for all w e Fc. 

Define characters G{'Hn) C "H^ and infinitesimal characters q{'Hn) C Wj^ as 

GCHn) = {a eUl,: a{I) = I, a{LJiUJu!2) ^ a{uji)a{uj2) ioT uJi,uj2 eFc} (48) 

giUN) = {aeUl,: a{I) = 0, a{ujiinLU2) = for LUi,uj2eFc\{l}}. (49) 

The convolution product on H^ is defined in the standard way: 

ao/3{uj) = ^ a(w(i))/3(cj(2)), (50) 

using the Sweedler notation. The convolution is the extension of the planar Grossman-Larson 
product from finite series to infinite series by considering H*j^ as the projective limit H*j^ ~ lim A/fc, 
where A4 = spanjo; g F: |a;| < k}. Note that the series are formal power series, and convergence 
in concrete cases, such as flows on manifolds, must be addressed by additional theory. 

GCHn) with the convolution product forms a group called the character group of Hn, where the 
unit and the inverse is given by the unit and the antipode in the Hopf algebra Hn, see [22) . In the 
special case where the post-Lie algebra is pre-Lie, this is the Butcher group, first introduced in [1] as 
a tool to study numerical integration. More generally, the elements in qIT-Ln) can represent vector 
fields and elements in G{Hn) diffeomorphisms on a manifold M. The convolution represents the 
composition of diffeomorphisms. Parallel transport of g S H^ along the t = 1 fiow of / G q(T-Ln) 



is represented by the exponential of the connection, which using ( 42 ) becomes 



exp(/[>)g ■= g + f '> g + -f > {f '> g) + ■ ■ ■ = {1 + I + -fof + •••)> 5 = cxp° ft>g, 
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where exp° is the exponential with respect to the composition product. The map cxp° : g('HAr) — > 
G{T-Ln) is 1-1, with the inverse given by the eulerian idempotent }22| . 

We win now be more concrete and discuss basic apphcation of Li3-series in numerical analysis 
of integ ration on a homogeneous space A^. Consider the post-Lie algebra on AA introduced in 



Section 2.2.3 written compactly as g^ :— il'jj(G,Q) ~ T{G Xh q), and define the enveloping 
algebra U{g)^ := fl%{G, U{q)) ~ T{G Xh U{g)). Let a: M ^ G he a. Cartan gauge (left partial 



moving frame), with a^ : TAi — > G Xh g and p* : G Xfj g — > TM as defined in Section 2.2.3 
Suppose we want to integrate a differential equation on M given as 

y'(t)=F(y(i)), 2/(0)= 2/0 eX, 
where F £ XM.. We represent the vector field as i^ = p*{f), where 

A numerical method with timestep < ft. G M is a diffeomorphism 'i'hf- M ^ A4, which is 
not identical to the exact flow ^hf (but preferably quite close). Many numerical methods are 
maps which can be represented as LB-series, and the analysis of the LB-series of the numerical 
solution is a fundamental tool to answer many questions about accuracy and geometric properties 
of the numerical algorithm. We will not go into details about this here, but merely discuss some 
different ways a LB-scries can be interpreted as a flow ^1^;^^ on A4. Consider "H^, where C — {•}, 
and an identification • iH- / e q^ . This extends uniquely to a map Tf-. H*^ -^ U{q)^. For a 
forest uj £ Fq, J^f{oj) is called an elementary differential operator. Note that for t g M we have 
IFtfii^) — ^'"'-^/(w). B-series and LB-series are traditionally considered as time-dependent series 
of differential operators on M" and on Ai respectively, for a G H*j^ given as 

wGFc 

There are (at least) three different ways LB-series are used to represent the numerical method. 

Parallel transport. Find a £ G{Hn) such that 

Thfia)t>g^nf9, for aU g e C/(G)^, (51) 

where ^Jj^g denotes parallel transport of g along "^hf- 

Backward error. Find f3 E g(HAr) such that '^hf is exactly the t = 1 flow of an autonomous 
vector field Fh — p^,fh, where ^ 

fh = ThfW)eg^. (52) 

Development. Given a curve t i— > y{t): / C M i— >■ A^ and a left (possibly partial) moving 
frame cr: A^ ^- G with Darboux derivative uja € il^{A4,Q). The curve S: I C M. —>■ q given as 
S{t) — uja-oy'{t) is called the development oiy(t) with respect to the moving framerl A computation 
yields: 

Proposition 4.2. Let y[t) E Ai be an integral curve of a differential equation 

y\t)=F{y{t)), y(0) = 2/o. (53) 

The development of y with respect to a is S{t) ~ f{a{y{t))), where f = a^,F G g^ ~ ^^u{G,q). 
Equivalently 

Sit) = ^IjfHyo)), (54) 

where ^*~/ denotes parallel transport of f along the flow of itself. 
'■J 



^Example: The classical Frenet-Serret frame of a curve in R'^ is the development with respect to the full moving 
frame obtained by prolongation of the Euclidean motion group acting on R^ . 
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We can represent the curve in y{t) in (53) algebraically by finding 7 £ 0('Hjv) such that 



m = ^-^t/(7)(yo) (55) 



is the development of y{t). 



Remark 4.3. The three ways of representing the numerical method, as a e G{Hn) representing 
parallel transport along the numerical flow ^/i/, as /3 € 0('HAf) representing the backward error 
vector field F and as 7 G q{'Hn) representing the development of the integral curve of the numerical 
flow, are all discussed in [55] . However, the third representation 7 is presented differently, and the 
close relationship between this representation and moving frames is therefore not obvious. 

The algebraic relationship between a, /3 and 7 is discussed in [22]. We find a from jj by 
applying exp° and /3 from a via the eulerian idempotent in 'Hn- From a we find 7 by applying 
the Dynkin idempotent in T-Lni and conversely a is found from 7 by a formula involving certain 
non-commutative Bell polynomials. 

Remark 4.4. Consider the exact flow ^tf of the differential equation y' = F{y). In this case, 
obviously, /? = •. From this we can compute explicitly the development S{t) G g of the solution 
curve y{t) G M, see [22], yielding 



m 





For instance, given a curve y{t) G M'^ where the moving frame is obtained by prolongation of the 
euclidean motion group acting on K^, this expresses the Taylor expansion of the classical euclidean 
curvature and torsion of the space curve y{t), in terms of elementary differentials of /. 

Concluding remarks 

The theory in this paper opens several interesting areas of further investigation. We are convinced 
that post-Lie algebraic structures will find applications in many areas also outside numerical 
analysis, such as stochastic differential equations, control theory and sub-Riemannian geometry. 
Within numerical analysis, the paper points to different ways of applying moving frame techniques 
in geometrically 'nice' ways of choosing isotropy in Lie group integration. This still has to be 
investigated numerically and computationally. 
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